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ABSTRACT 



We report on the results of two coordinated multiwavelength campaigns that focused 
on the blazar Markarian 421 during its 2006 and 2008 outbursts. These campaigns ob- 
tained UV and X-ray data using the XMM-Newton satellite, while the gamma-ray data 
were obtained utilizing three imaging atmospheric Cerenkov telescopes, the Whipple 
10m telescope and VERITAS, both based in Arizona, as well as the MAGIC telescope, 
based on La Palma in the Canary Islands. The coordinated effort between the gamma- 
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ray groups allowed for truly simultaneous data in UV/X-r ay /gamma-ray wavelengths 
during a significant portion of the XMM-Newton observations. This simultaneous cov- 
erage allowed for a reliable search for correlations between UV, X-ray and gamma-ray 
variability over the course of the observations. Investigations of spectral hysteresis and 
modeling of the spectral energy distributions are also presented. 

Subject headings: BL Lacertae objects: individual (Markarian 421) — galaxies: active 
— gamma rays: observations — radiation mechanisms: non-thermal — X-rays: galaxies 



Introduction 



Blazars, a sub-class of active galactic nuclei (AGN), are some of the most intriguing sources 
in the high-energy sky. Their rapid variability and the non-thermal nature of their emission, 
presenting a continuum across nearly the entire electromagnetic spectrum, imply that the observed 
photons originated with in highly relativistic jets oriented very close to the observer's line of sight 
(jUrrv Padovanil Il995l ). Therefore, blazars are excellent laboratories for studying the physical 
processes within the jets of AGN. They were among the first very high energy (VHE; E >100 GeV) 
sources to be detected and today there are 25 known VHE blazars, including flat spectrum radio 
quasars and BL Lac objects. 

Various models have been proposed to account for the broad-band spectral energy distributions 
(SED)s observed in VHE blazars, which typically display a characteristic double peak when plotted 
as vF v against u, with peaks occurring at keV and TeV energies. The models are generally divided 
into two classes: leptonic and hadronic. Both leptonic and hadronic models attribute the peak at 
keV energies to synchrotron radiation from relativistic electrons (and positrons) within the jet, but 
they differ on the origin of the TeV peak. The leptonic models advocate the inverse Compton scat- 
tering mechanism, utilizing synchrotron self Compton (SSC) interactions and /or inverse Compton 
interactions with an external photon field, to explain the VHE emission (e.g., iMarscher Gear 
19851 : iMaraschi et al.lll992l : iDermer et al.lll992l : ISikora et al.lll994l ). The hadronic models, however, 
account for the VHE emission by 7r° or charged pion decay with subsequent synchrotron and/or 
Com pton emission from decay products, or synchrotron radiatio n from ultra-relativistic protons 
(e.g., lMannheimlll99.4 lAharonianl[200nl : IPohl fc Schlickeiserlhood ). 



Observationally, blazars are known to undergo both major outbursts on long time scales and 
rapid flares on short time scales, most prominently at keV and TeV energies. During some outbursts, 
both of the S EP peaks have been obs erved to shift towards higher energies in a generally correlated 
manner (e.g.. iBlazejowski et al.ll2005l ). The correlation of the variabilities at keV and TeV energies 
(or lack thereof) during such outbursts has aided in refining the emission models. In addition, rapid, 
sub-hour flaring is interesting as it provides direct constraints on the size of the emission region. 
These rapid flares also present an observational challenge to multiwavelength studies, as truly 
simultaneous data must be used in order to develop a reliable characterization of the broadband 
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behavior of these objects. 



Markarian 421 (Mrk421; 1101+384), at a redshift of z = 0. 031, was the first blazar, as well 
as the first extragalactic source, to be detected at TeV energies (jPunch et al.lll992l ) and has since 
remained one of the most active VHE blazars. Its SED has peaks at keV and TeV energies, and 
it has been known to dem onstrate rapid, sub-hour, fla ring behavior at these energies during the 
course of an outburst (e.g.. lCuill2004! ; lGaidos et al.lll99q ) indicating very compact emission regions. 
Many extensive multiwavelength campaigns studying Markarian 421 have been undertaken at keV 
and TeV energies during outbursts, but the degree of simultaneity of the multiwavelength coverage 
varies and is often not adequate to account for the most rapid variability of the source. Due to 
these considerations, one must exercise caution in the interpretation of some of the results from 
these campaigns. 

In this work, we report on results from a target-of-opportunity (ToO) program on VHE blazars 
that makes use of a unique combination of capabilities provided by the XMM-Newton satellite and 
several ground-based imaging atmospheric Cerenkov Telescopes (IACT)s to obtain truly simulta- 
neous coverage in the optical/UV, X-ray and VHE bands. Unlike other satellites frequently used 
for multiwavelength campaigns, XMM-Newton 's highly elliptical orbit allows for long observations 
that are not frequently interrupted by earth occultation. In addition to its X-ray instruments, 
XMM-Newton carries an optical/UV telescope co-aligned with the X-ray telescopes and can thus 
provide simultaneous coverage in the optical/UV band. While our primary emphasis is on the 
X-ray and VHE bands, the optical/UV coverage helps to constrain the overall SED shape. 



2. Observations and Data Reduction 

Our ToO program was first triggered in April 2006 by a major outburst from Markarian 421 
as detected by regular monitoring of the VHE band by the Whipple 10m telescope. Because of 
XMM-Newton visibility constraints, the coordinated multiwavelength observations did not take 
place until after the peak of the outburst as indicated by the overall monitoring campaign. The 
program was triggered again in May 2008 by another major outburst from Markarian 421 detected 
in the VHE band. However, we once again captured only the decaying portion of the outburst. 
However, taken together, the two campaigns have produced a significant amount of simultaneous 
optical/UV, X-ray and VHE data on the source. 



2.1. XMM-Newton Observations 



The X-ray and optical/UV obse rvations were taken by the XMM-Newton satellite's EPIC-pn 
(EPN) detector (jStriider et al.l 120011). covering a spectral range of approximately 0.5 - 10 keV, and 



the Optical Monitor (OM; iMason et al. Il200ll ). capable of covering the range between 170 and 650 



nm (7.3 eV and 1.9 eV). The metal oxide semiconductor (MOS) detectors were also operated during 
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some of the observations, but the data were not used here. During both observations, the EPN was 
operated in fast timing mode to minimize photon pileup at the expense of imaging capability along 
the direction of event readout. The thin optical filter was in place for both observations. The 2006 
and 2008 observations produced EPN exposures of approximately 42 ks and 43 ks, respectively. 
Due to the brightness of the source, during the first half of the 2008 observations both the EPN and 
MOS detectors entered into counting mode, resulting in a loss of frames. In an attempt to avoid 
further frame losses, the MOS detectors were shut down to provide more bandwidth for telemetry 
to the EPN detector. As a result of the detector being in counting mode, the first half of the EPN 
observation suffers from telemetry gaps which can be accounted for and corrected as explained 
below. 

In parallel with the EPN observations, a series of exposures was taken utilizing the OM. The 
2006 and 2008 observations produced 15 and 10 exposures, respectively, with total exposure times 
of 32.5 ks and 22.0 ks. For all exposures, the OM was in imaging mode with the UVM2 filter (200 
- 300 nm) in place. Table [1] summarizes these observations. 



The 2006 EPN data were initially processed using XMMSAS v7.0 (jGabriel et all 12004 ). Stan- 
dard XMM-Newton data analysis procedures were followed to filter and reduce the data, generate 
the various lightcurves, extract the source and background spectrum, and generate the RMF and 
ARF files for subsequent spectral analyses. In summary, the event list was examined to check 
for periods of soft proton flaring, where background event rates are higher than 0.4 cts/sec (for 
events with energies E > 10 keV). In addition, "bad events" were removed including events which 
were close to CCD gaps or bad pixels, and only single and double photon events were included 
(PATTERN<4). A one dimensional histogram of counts was produced to determine the position 
of the source. Based on the histogram, we chose columns 30-45 (in RAWX) for the source region 
and columns 1-15 and 57-64 for the background region. These regions were used to extract source 
and background spectra. 

The 2008 EPN data were processed similarly, using the same filtering criteria, but with XMM- 
SAS v8.0. In this case the source region was determined to be in columns 31-45 and the background 
regions were taken to be columns 1-19 and 57-63. It should be noted that a new SAS task, epiclc- 
corr, was used to correct for the effects of deadtime, telemetry saturation, and frame dropouts in 
producing the light curves from the 2006 and 2008 data. 

The OM data were analyzed with the standard SAS task omichain, which produces images 
and source lists for all OM exposures present in any given observation. The 2006 OM observations 
presented here are a combination of 15 exposures with 2000 s and 2500 s integration times while 
the 2008 data comes from 10 exposures each with a duration of 2200 s. 
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2.2. VHE Gamma-Ray Observations 

Three ground-based IACT facilities were utilized to observe the source in the VHE band 
during the XMM-Newton observations. These IACTs detect gamma-rays by imaging the flashes 
of Cerenkov light emitted by gamma-ray induced electromagnetic showers within the atmosphere. 
The Major Atmospheric Gamma Imaging Cerenkov (MAGIC) telescope provided 14 ks of coverage 
during the first part of the 2006 XMM-Newton observation, and the Whipple 10m telescope provided 
12 ks of coverage during the latter part of the 2006 observation. The geographic separation of the 
experiments allowed for extended, simultaneous coverage during the XMM-Newton observation. 
During the 2008 XMM-Newton observation, only the Very Energetic Radiation Imaging Telescope 
Array System (VERITAS) was used, providing 9 ks of VHE coverage. 



2.2.1. Whipple Data 



The Whipple 10m IACT (jKildea et al.l 120071 ) focuses Cerenkov light on a camera composed 



of 499 photo multiplier tubes (PMTs). The energy threshold for the Whipple 10m is near 400 
GeV. There are two modes used for taking data. During the 2006 observation, two hours of data 
were taken in TRACKING mode, where the source is centered in the camera, and two hours of 
data were taken in ON/OFF mode, where half of the data runs are offset from the source by 30 
minutes in right ascension to provide an independent measurement of background events (primarily 
cosmic-rays). The runs each had a duration of 28 minutes, and the total time on source for the 
Whipple 10m telescope was 3 hours on 2006 April 30 from 4:37 UT to 8:05 UT. The observations 
were taken with a source elevation range of 45° - 79° , with an average elevation of 64° and during 
favorable weather conditions. 

Whipple data reduction involves two stages, image cleaning and image parameterization. In 
the cleaning stage, the data is flat-fielded using a run during which a nitrogen arc lamp is pulsed 
to illuminate the camera's pixels. Any differences between night-sky background between ON and 
OFF runs is also accounted for at this stage. In addition, pixels that make up shower images 
are s elected and the rem aining background pixels are then removed prior to image parameteriza- 
tion (jPunch et al.l Il99lh . During the image parameterization stage, each Cerenkov image in the 



telescope's camera is characterized by an ellipse, using a moment analysis of the recorded signal am- 
plitudes in each pixel. Each shower is characterized by the major axis {length), minor axis {width), 
the angle between the major axis of an image's ellipse and a line from the centroid of the ellipse 
to the source position (a), the length of this line {distance), and the overall signal of the shower 
{size). These are known as Hillas parameters and they allow for the removal of the primary source 
of background, cosmic rays, by expl oiting the in trinsic differences in the development of hadronic 



cosmic ray and gamma-ray showers (lHiilaslll985l ). These differences result in different distributions 
of light at ground level and cuts can be used to select events most likely to be gamma-rays. Stan- 
dard cuts were made on the Hillas parameters, as shown in Table El to reject as many background 
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events as possible while still retaining many gamma-ray candidates. 

The runs taken in ON/OFF mode obtain a direct measurement of background from the off- 
source runs, while in TRACKING mode a "tracking ratio" is used to estimate the background. 
Gamma-ray events should have a small a parameter if the source is centered in the field of view, 
so events with large values for a can be used to estimate the background if the ratio between 
background rates at small and large values of a is known in the absence of a source. This "tracking 
ratio" used here was found by taking observations of blank sky fields with no detected VHE sources 
throughout the 2005 - 2006 observing season. It is defined as the ratio of the integrated number 
of events between a = - 15° and those between a = 20 - 65°. Once it is calculated, the tracking 
ratio is used to estimate background rates for all TRACKING observations. 



2.2.2. MAGIC Data 



MAGIC is the largest single dish Cerenkov telescope in operation (jBaixeras et al.ll2004l : ICortina et al 



2005). A 17 m tesselated reflector focuses the light from air showers on a camera composed of 576 
PMTs. For high elevation angle obs ervations, the MAGIC trigger threshold currently reaches down 
to 50 - 60 GeV d Albert et allhoosh . 



The measurements reported in this article were conducted from 2006 April 29 21:32 UT to 
2006 April 30 00:59 UT at elevation angles spanning 4 9° - 81°. MAGIC observed the source 
employing the so-called wobble mode (jFomin et al.lll994l ). during which the telescope alternates 
between tracking two (or more) opposing sky directions, each 0.4° off-source, for 20 minutes each 
(jAnderhub et al.l |2009j) . After removing events not containing sufficient information for further 
analysis (see Table [3] top), events from accidental triggers, triggers from nearby muons, and data 
affected by adverse meteorological conditions, 2.8 hours out of the 3.8 hours of data were used for 
further analysis. 

The data were process ed using the analysis and reconstruction software package for MAGIC 
data dBretz fc Domeil2008h . A d escript i on of t he different analysis steps can be found in lGaug et al.l 



data ( B 



Albert et al 



(2008J) and iBreta ([2005). For this analysis, the signal has been extracted 
using a spline algorithm. After calibration, the shower images were cleaned of background noise by 
requiring a minimum photoelectr on signal in the pixels as well as temporal coincidence with adjacent 
pixels (time image cleaning, see lAnderhub et al.l (|2009i )). The recorded events are characterized 
by several image parameters based on the shower light distribution (amongst them the previously 
me ntioned Hillas paramet ers) and on the temporal shower development in the camera plane {slope, 



see lAnderhub et al.l (120091 ) ). True gamma-ray events coming from the observed source are extracted 
from the hadronic background by cuts in this image parameter space, using the separation cuts 
listed at the bottom of Table El Events originating from the source are selected by a cut on 
9 2 (where 9 is the angular distance between the expected source position and the reconstructed 
gamma-ray arrival direction). Gamma-ray events are then separated from background events by a 
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size dependent parabolic cut in width x length x n (jRiegel et alll2005l ). 

The Markarian 421 observations presented here are among the first data taken by MAGIC after 
major hardware updates in April 2006. Due to these changes, a thorough examination of the data 
was undertaken to ensure its reliability. Despite the hardware changes, all MAGIC subsystems were 
working as expected with the exception of an unstable trigger behavior for some PMTs, leading to 
an early signal arrival time for those affected. A careful study of the systematics in the low energy 
region has been condu cted, leading to the decision of raising the energy threshold to 250 GeV 
(jAnderhub et alibood ). 



2.2.3. VERITAS Data 

VERITAS is an array of four 12 m diameter IACTs each focusing light on a camera consisting of 
499 PMTs. Utilizing multiple telescopes provides a stereoscopic view of showers, making it possible 
to significantly reduce muon events. This substantially improves the low energy performance of 
the detector since images produced from nearby muons are otherwise virtually indistinguishable 
from images produced by low-energy gamma-rays. Therefore the array allows for a lower energy 
threshold than possible with a single telescope of the same size operating alone. For VERITAS 
this threshold is near 275 GeV. In addition, stereoscopic observations facilitate a much improved 
determination of the core position of showers and thus provide improved energy resolution and 
background rejection. 

During the 2008 observation, VERITAS took data on Markarian 421 while operating in wobble 
mode, with an offset angle of 0.5°. To provide background estimatio ns using this observation mode, 



we follow the reflected-region model (jBerge. Funk. &: Hintonl 120071 ) where one collection region is 



placed at the source position and others of equal size are placed at equal offsets from the center of 
the field of view and used for background measurements. VERITA S observed Markarian 421 for 2.5 
hours on 2008 May 07 from 3:59 UT to 6:28 UT using a series of 20 minute runs, with an effective 
exposure time of 129.7 minutes. The observations were taken as the source elevation ranged from 
79° to 55°. Data quality checks confirmed good weather during the observations. 

The VERITAS data from 2008 was analyz ed with Veritas Ga mma-ray Analysis Suite (VE- 



GAS), the standard VERITAS analysis package (jCogan et al.ll2007l ). The shower images were first 



corrected for relative gains and cleaned to remove isolated pixels. The images were then parame- 
terized in a similar manner as for Whipple analysis. Initial quality cuts were made on the images' 
size, pixel count and distance. For this analysis, events that survived the initial cuts but only 
contained data from the two telescopes with the smallest physical separation were removed. This 
smaller separation distance results in a greater number of muons surviving the initial cuts, and 
in addition these events can distort the calculation of impact distance due to the relatively small 
distance between the telescopes. After stereoscopic reconstruction and calculation of each event's 
impact distance (distance from the shower core to the telescope), the events were parameterized 
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using the quantiti es mean scaled width (MSW) and mean scaled length (MSL) (|Konopelkolll995l ; 



Daum et al.lll997l ). These parameters are found by scaling the width and length parameters for 



a given telescope by the average expected values from simulations given an impact distance and 
size and then finding the average for all telescopes involved in an event. Cuts on MSW and MSL 
were made to separate gamma-ray events from cosmic-ray events. Table [3] details the cuts that are 
discussed above. A summary of all VHE gamma-ray data is shown in Table [5j 



Results 



3.1. Time Averaged VHE Gamma- Ray Spectra 



We followed the procedures from lMohantv et al.l (|1998l ) in constructing a spectrum based on 
the Whipple data. This method depends on separate energy estimates for both on-source and off- 
source runs, so the runs taken in the TRACKING mode were matched to contemporaneous OFF 
runs taken at similar elevation angles. The selected matched runs were taken within two days of the 
Markarian 421 observations and at elevation angles within 1° of the corresponding TRACKING 
runs. The energy spectrum was fit using a simple power law: 



dN 
~dE 



F ■ 10 



-ii 



— v 

ITeV/ 



photons 
TeV cm 2 s 



(1) 



Finding best-fit parameters of a = 2.23 ± 0.38 and Fq 
(P = 98.2%). Uncertainties are statistical only. 



2.65 ± 0.77 yielding a x*/d.o.f. of 0.17/3 



The MAGIC spectrum has been derived on the same data basis as the light curve but a 
somewhat looser area cut was applied. This cut yielded a constant cut efficiency as a function 
of energy of 90% for Mo nte Carlo simulate d gamma-ray events, increasing the gamma-ray event 
statistics at the threshold ([Albert et al.ll2009i) . The energies of the gamma even ts were reconstructed 
using a random forest regression method (|Breimanl l200ll ; lAlbert et al.l l2007d ) trained with Monte 
Carlo events. The MAGIC spectrum was modeled similarly to the Whipple spectrum, finding best 
fit parameters of a = 2.28 ± 0.09 and F = 2.49 ± 0.17, yielding a x 2 / d -°- f - of 2.04/4 (P = 72.9%). 
All the stated uncertainties for MAGIC are purely statistical. The energy scale is known with an 
uncertainty of ±16%, the flux normalization within a systematic error of 11% (not including the 



energ y scale error), and the fitted power law slope has a systematic uncertainty of ±0.2 ([Albert et al 
20081 ). 



For the MAGIC analysis, an additional systematic uncertainty is provided to account for 
possible effects arising from the hardware instability mentioned above. These effects consist of 
a moderate loss of low-energy showers as well as a minor additional uncertainty in the image 
parameter calculation for showers of higher energy. The effect on the differential flux level is 
estimated to be 10% from 250 GeV to 400 GeV and 3% for higher energies (jAnderhub et al.ll2009i ). 
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The spectral a nalysis for the VERITAS data was performed using the VEGAS analysis package 
(jCogan et al.ll2007l ). with the same model for background estimates as used for the light curve. The 
spectrum was again fit with a simple power law model, finding best fit parameters of a = 2.91 ±0.13 
and Fq = 2.01 ± 0.15, yielding a x 2 /d.o.f. of 9.66/8 (P = 29.0%). Again, stated uncertainties are 
statistical only. The results from the spectral fits for all VHE data can be seen in Figure [TJ 



3.2. Spectral Energy Distribution and Modeling 

The X-ray spectrum was initially fit using XSPEC 12. The data were fit with a power law 
modified by interstellar absorption, yielding a value for the photon index of a = 2.258 ± 0.002 
and 2.153 ± 0.002 for the data obtained simultaneously with the 2006 MAGIC and Whipple 10m 
observations, respectively. For X-ray data taken during the 2008 VERITAS observations, a photon 
index of a = 2.519 ±0.010 was found. The hydrogen column density was left as a free parameter for 
all fits, finding values of 3.16 x 10 20 cm~ 2 , 2.25 x 10 20 cm~ 2 and 4.51 x 10 20 cm~ 2 for the data obtained 
simultaneously with the MAGIC, Whipple 10m and VERITAS observations, respectively. Using 
these results, the spectrum was unfolded and de-absorbed to derive the intrinsic X-ray spectrum of 
Markarian 421. Only statistical errors were taken into account here. In addition, the count rates 
found for the OM exposures were converted to flux using the standard conversion factoi0 and an 
average p oint was determined for each time interval. In addition, using the ultraviolet extinction 



law from ICardelli et al.l (|1989l ). the absolute extinction for the UVM2 band was calculated to be 



A(UVM2) = 0.13, allowing for de-reddening of the OM data using a correction factor of 1.13. 

Figure [2] shows the broadband SEDs corresponding to the three epochs of VHE observations. It 
is important to note that within each epoch the multiwavelength data are genuinely simultaneous. 
Spectral variability is observed between epochs. 



Modeling of the SEDs was carried out using a leptonic model (jBottcher Chian g 2002). In 
this model, the spectral distribution of injected electrons is described by a power law with low 
and high energy cutoffs of 7 m j n and ~f ma x, respectively. The emitting region is assumed to be in 
a state of temporary equilibrium and of spherical shape, with radius R, and moves out along the 
jet at relativistic speed v/c = (1 — 1/T 2 ) 1//2 , where T is the bulk Lorentz factor. As the emitting 
region moves along the jet, particles cool due to radiative losses and may escape from the region. 
The timescale of these escapes is factored into the model as t esc = rjR/c, with 77 > 1. The radia- 
tive processes considered include synchrotron radiation, SSC, and inverse Compton scattering of 
external photons. However, we found that a model with a negligible contribution from external 
photons (i.e., a pure SSC model) provides a good match to the SEDs during both observations. 
In addition, the SEP matches have b een absorbed with the extragalactic background light model 



discussed in iFranceschini et al.1 ((2008). The values for magnetic field and the spectral index of the 



See: 



http://heasarc.nasa.gov/docs/xmm/sas/USG/nodel35.html 
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injected electrons were varied until a good match was found for the 2006 and 2008 SEDs. The 
parameters for these models are shown in Table [6) The models are also shown in Figure [2j 



3.3. Cross-Band Correlation 

The VHE, X-ray and UV light curves are shown in Figure [3] with 1-a uncertainties. The X-ray 
lightcurve was initially binned using a time interval of 500 seconds, while the bins for the VHE light 
curves were primarily determined by the standard length of data runs for the respective telescopes. 
The Whipple 10m, MAGIC, and VERITAS data are in 28-minute, 22-minute and 20-minute bins, 
respectively. To provide a direct comparison between the results obtained with different VHE 
experiments, we show the gamma-ray fluxes above a common energy threshold (250 GeV). To 
reach this threshold for the Whipple 10m, the flux was extrapolated from the power law spectrum 
fit to the data. 

In addition, contemporaneous data taken on the Crab Nebula, a standard candle for VHE 
experiments, were studied to assess systematic uncertainty in the flux calibration of the VHE data. 
The systematic uncertainty in flux calibration between the MAGIC and VERITAS experiments 
was estimated to be on the order of 10% with MAGIC systematically measuring a lower flux than 
VERITAS. The systematic uncertainty in the Whipple 10m and VERITAS flux calibration was 
found to be on the order of 40% with the Whipple 10m systematically measuring a higher flux than 
VERITAS. These systematic uncertainties are not included in the data. 

The error bars on the VHE data are relatively large, so in order to quantify the variability, 
we first tested the entire light curve by fitting a constant flux value to the data. The resulting 
fit had a xi °f 3.92 for 20 degrees of freedom. Similarly, as the MAGIC data made up the most 
variable portion of the VHE light curve, these data alone were fit to a constant, with a Xu °f 5.22 
for 6 degrees of freedom. This inconsistency of the data with fits to a constant indicates significant 
variability in the VHE band. Markarian 421 also varied significantly at X-ray energies during 
the 2006 observation, with the count rate initially decreasing during the course of the MAGIC 
observation and slowly increasing during the Whipple observation. Though both showed significant 
variability, the X-ray and VHE data from 2006 do not appear to be correlated. 

To examine the X-ray /VHE correlation more closely, we show in Figured the measured VHE 
flux and X-ray count rates. The X-ray data were rebinned to match the resolution of the corre- 
sponding VHE data. To quantify the correla tion and examine th e uncertainty in this correlation, 



we used a method similar to that discussed in lAlbert et al.1 (|2007bl ). A set of 25000 light curves was 
simulated based on the Gaussian errors of the data points in the X-ray and VHE bands, respec- 
tively. For each pair of simulated X-ray/VHE light curves, the value of Pearson's r was calculated. 
A histogram of these possible r values was generated, giving an average r value of —0.050 ± 0.050, 
indicating a lack of correlation. 
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To compare to previou s work, the best fit correlation between the X-ray and VHE bands found 
tazeiowski et al.l (|2005l ) is plotted in Figure H (dotted line). In order to carefully compare this 
previous result to our results, the PIMMS toojj was used to convert the RXTE count rates to 
XMM-Newton count rates, taking into account energy range, spe c tral s hape and hydrogen column 
density. To more directly compare our data to iBlazeiowski et al.l (|2005l ) , we scaled the correlation 
to match our Whipple 10m points in Figured] (dot-dashed line). The overall normalization shift 
could be a reflection of hysteresis on long timescales. Although the scatter in the VHE points are 
large, the VERITAS points are systematically below the scaled correlation. This effect is made 
worse by the fact that the Whipple 10m points are known to be systematically higher than the 
VERITAS points (by 40%). 

There appears to be no correlation between the VHE and UV variations, most notably during 
the first half of the 2006 observation, where the UV rates increase as the VHE flux decreases. 
However, the UV rates during both observations appear to roughly follow the trend of the X-ray 
rates, with a significantly higher rate during the 2008 observation. Using the same method as used 
above to examine correlations between the VHE and X-ray bands, we found an average value for 
Pearson's r of 0.940 ± 0.001 for the total set of X-ray and UV data, indicating strong correlation. 



3.4. Spectral Hysteresis 



We also examined the spectral evolution of the X-ray data to determine if any hysteresis was 
present during the observations indicating a dependence of the system on previous states. The data 
were divided into three energy bands in order to calculate hardness ratios: 0.5 - 1 keV, 1-3 keV, 
and 3-10 keV. These bands were chosen such that each band contains a roughly equal number of 
counts. Figure [5] shows the evolution of the hardness ratios versus intensity (in terms of count rates) 
through the flare seen in the 2006 data (which peaks near 17-18 ks into the observation; see Figure 
[3]). Spectral hysteresis is clearly seen during this rapid flare. Clockwise patterns indicate a lag in the 
response of lower-energy photon s with respect to that o f higher-energy photons. Similar patterns 
have been seen previo usly (e.g., iBrinkmann et al.ll200ll ). though counter-clockwise patterns have 
been observed as well ()Cuill2004j ). Different hysteresis patterns indicate a complex mechanism that 
may differ between outbursts. 

No rapid flares were observed during the 2008 observation. Examination of the entire 2008 
observation did not reveal any overall spectral hysteresis pattern occurring on the order of several 
hours. Similarly, no overall hysteresis was found when considering the entire 2006 observation. 



2 See: |http://heas arc.gsfc.nasa.gov/ Tools/w3pimms. html 
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Discussion 



This work provides the results from a ToO campaign with the primary focus of studying the 
rapid flaring activity of blazars on sub-hour timescales in both the X-ray and VHE bands. In 
addition, these observations expand the pool of truly simultaneous, multiwavelength observations 
of blazars, essential for the study of correlations during rapid flares. The variability observed during 
the 2006 outburst in the simultaneous X-ray and VHE observations should be enough to provide 
useful information about a correlation, if any, between the two energy bands. In addition, the 
observations from 2008 during a separate outburst, provide an extended picture for investigating 
possible correlations between these two bands. Rapid X-ray and VHE flares are often expected to 
be correlated in one-zone SSC scenarios, as all photons are expected to originate from the same 
population of electrons. Examination of the observed SEDs can provide additional constraints on 
model parameters and reveal what factors strongly affect each flare's spectral profile. 

The SED fits show a steeper spectrum for the 2008 observation than for the 2006 observation, 
particularly at X-ray energies. This is intriguing; as previous investigations have observed spectral 
hardening with increased flux (e.g., Kue et al.1 120061 ). In addition, both peaks in the SED show a 
slight shift to lower energies during the 2008 observation where previous investi gation has shown 
a shift to higher energies with increased luminosity (e.g., IBlazejowski et al.ll2005l ). Figure [2] shows 
that our one-zone SSC model fits both the 2006 and 2008 data quite well. 

Historically, multiwavelength monitoring campaigns have observed correlation between the X- 
ray and VHE emission in Markarian 421, though the data sets are not comprised enti rely of strictly 
simultaneous multiwavelength data (e.g.. IBlazejowski et al.ll2005l : iFossati et al.ll2008l ). Though the 
dynamical range of the observed variability is not as large as that seen in previous long-term 
monitoring campaigns, we seem to be seeing short-term variability that shows a different correlation 
pattern here. Surprisingly, both the 2006 data alone and the entire 2006/2008 data set show no 
obvious correlation, with a Pearson's r for the entire X-ray /VHE data set of —0.050 ±0.050, though 
the SSC model fits the SEDs well (see Figure[2]and Figured]). This implies that the X-ray and VHE 
photons may originate fr om electrons with diffe rent energies, similar to what was observed during 
a study of PKS2155-304 (jAharonian et al.l 12003 ) . Other scenarios that could explain the observed 
variability patterns include the possibility of an inhomogeneous emission region or hadronic origin 
of the VHE emission. It is interesting to note that Marka rian 421 is clearly behav ing very differently 
here from what is usually reported of the source (e.g., IBlazejowski et al.ll2005l ). where the X-ray 
and VHE variabilities are seen to be strongly correlated. Such personality makes it difficult to 
generalize the results to other blazars. 

In addition, the OM allowed for a detailed search for UV/X-ray or UV/VHE correlations. 
Though previous multiwavelength campaigns on Markarian 421 have obtained optical, X-ray and 
VHE data, correlati ons between the X-ray and optical data have either not been st udied or were 



not significant (e.g.. lAlbert et al.ll2007al : IBlazejowski et al 



2005 



Horan et al.l 12003 ). This study 



found the average rate from the XMM-Newton PN detector during the 2008 observation was about 
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20% higher than the rates measured during the 2006 observation. Examination of the OM data in 
the UV band shows that the count rate more than doubles between the 2006 and 2008 observation. 
The UV rate also appears to vary in step with the X-ray rate within the individual observations. 
Calculation of the Pearson's r value for the total UV/X-ray data set yielded a value of 0.940±0.004, 
indicating a strong correlation. This provides direct observational evidence for a link between the 
emission mechanisms at X-ray and UV wavelengths. On the other hand, we found no apparent 
UV/VHE or X-ray /VHE correlation, implying that the seed photons for the inverse Compton 
process (to produce the VHE photons) cannot be provided by the UV or X-ray emission observed. 

We also searched for hysteresis patterns in the X-ray data to provide further information 
about the outbursts. Hysteresis was observed during one rapid flare in the X-ray data from the 
2006 observation. No rapid flares occurred during the 2008 observation. Spectral hysteresis has 
been commonly observed i n blazars, but the phenomenon is not yet understood completely. In 
Kirk Mastichiadis (jl999l ). a simple model is used to produce spectral hysteresis patterns that 



may be observed during the course of a flare. In this model, the behavior is characterized by the 
relationship between three timescales associated with the duration of the flare variability (t var ), syn- 
chrotron cooling (tcooi) and particle acceleration (t acc ). The relationship between these timescales 
result in four possible cases that are discussed in detail. The clockwise hysteresis found in the X-ray 
data, indicating a lag at lower energies in the X-ray band, coupled with the essentially symmetric 
shape of the flare in the X-ray light curve seems to indicate that the case with t coo \ 3> t var t acc 
is most relevant to this observation. The lag in the low energy photons is a result of the inverse 
relationship between t coo i and the energy of the cooling particles. Although we observed clockwise 
patters in one of the observations, other patterns have also been observed. This is yet another 
example of the personality of the source and may indicate physical differences between individual 
flares and outbursts. 
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Fig. 1. — VHE spectral analysis results, showing the data obtained with MAGIC (squares, a = 
2.28 ± 0.09 and F = 2.49 ± 0.17), Whipple (diamonds, a = 2.23 ± 0.38 and F = 2.65 ± 0.77 ) and 
VERITAS (crosses, a = 2.91 ± 0.13 and F = 2.01 ± 0.15). 
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Fig. 2. — Spectral energy distribution with the SSC model for the 2006 (teal blue) and 2008 (red) 
data. The data from XMM- Newton OM, X MM- Newton EPN, MAGIC, Whipple and VERITAS are 
shown with triangles, filled circles, squares, diamonds, and crosses, respectively, with data taken 
during the MAGIC, Whipple and VERITA S observation times shown in green, blue and red. 
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Fig. 3. — Plot of light curves from (a) 2006 observations and (b) 2008 observations. XMM-Newton 
OM data (200-300nm) is shown with triangles in units of cts s _1 . XMM-Newton EPN data (0.5 - 
10.0 keV) is shown with filled circles in units of cts s _1 . Error bars for the EPN data are smaller 
than the data points. The data from MAGIC, Whipple and VERITAS are shown with squares, 
diamonds, and crosses, respectively in units of photons cm -1 s~ 2 (above 250 GeV). Note that the 
EPN and OM scales differ on the 2006 and 2008 panels. 
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Fig. 4. — Plot of simultaneous XMM-Newton EPN (0.5 - 10 keV) and VHE data to search for 
correlation between the two bands. MAGIC data points are designated by squares, Whipple points 
by diamonds, and VERITAS points by crosses. All XMM-Newton data used for this plot was 
analyzed with XMMSAS version 8.0.0 to m a intain consistency. The dotted line is the X-ray/ VHE 
correlation best fit from iBlazeiowski et al.l (|2005l ) converted to the appropriate units. The dot- 
dashed line has the same slope but is scaled to the average of our Whipple 10m results. 
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Fig. 5. — Hardness ratio vs intensity plots to study spectral hysteresis during the increase in 
activity centered around 17 ks into the 2006 XMM-Newton observation. Three bands are used for 
this analysis. The upper panel uses (3-10 keV)/(l-3 keV) for the hardness ratio, the middle panel 
uses (1-3 keV)/(0.5-l keV), and the bottom panel uses (3-10 keV)/(0.5-l keV). The arrows indicate 
the progression of time. Clockwise hysteresis is present in all three panels, indicating a lag in the 
lower energy bands. 
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Table 1: Summary of Markarian 421 XMM-Newton Observations. 

Obs. ID Start Time PN Exposure Avg. PN Rate Avg. OM Rate 

(ks) (cts s^ 1 ) (cts s" 1 ) 

0302180101 2006 Apr 29 4L9 336.30 ± 0.09 13.5±0.2 
20:44UT 

0502030101 2008 May 07 43.2 411.31 ± 0.12 35.82 ± 0.03 
01:34UT 



Table 2: Whipple Event Selection Parameters 



Parameter 


Values for Cuts 


Trigger level 


> 30 digital counts 


Shape cuts 


0.05°< width < 0.12° 




0.13°< length < 0.25° 


Muon cut 


length/size < 0.0004°d C - 1 


Distance cut 


0.4° < distance < 1.0° 


Orientation cut 


a < 15° 



Table 3: MAGIC Cuts. 



Parameter 


Values for Cuts 


Number of shower islands 


< 3 


Number of used pixels 


> 5 


Leakage 


< 0.3 


Gamma/Hadron Separation Cuts: 




e 2 


< 0.046 


Area 


< 0.265 x (1 - 0.0803 x (logW(size) - 5.77) 2 ) 


Standard time cuts 


slope > (distance — 0.5) x 7.2 



Note. - 



Leakage parameter is defined as the ratio of the signal in the outer pixels of the camera to the total signal. 
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Table 4: VERITAS Cuts 



Parameter 


Values for Cuts 


Size 


> 400 digital counts 


Pixels 


> 5 


Distance 


< 1.43° 


Mean Scaled Width 


0.05 < MSW < 1.16 


Mean Scaled Length 


0.05 < MSL < 1.36 
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Table 5: Gamma-ray data run information. 
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n /i no i nnn 
0.40 +000 


A A 

4.4 






TATA T 


OO 7 
ZO. / 


U T U I DDT T7 








onnr- A on 

2006 Apr 3D 


n A • ) — ' T TT 

04:37 U 1 


TT) A /^1T^T"\T/^1 


/i n 

4.9 




n r n r T tt 

05:05 UT 


1 RACKING 


5.4 




nr on ttt 

05:39 U 1 


1 RACKING 


6.6 




nr> nT TTTi 

06:07 UT 


1 RACKING 


3.7 




UD.O 1 U 1 








07:37 UT 


ON/OFF 


5.3 






TOTAL 


12.5 


VERITAS 








2008 May 07 


03:59 UT 


0.5° N 


15.5 




04:21 UT 


0.5° S 


15.4 




04:43 UT 


0.5° E 


16.9 




05:04 UT 


0.5° W 


13.5 




05:26 UT 


0.5° N 


15.2 




05:47 UT 


0.5° S 


14.6 




06:08 UT 


0.5° E 


13.7 






TOTAL 


39.6 



Note. - 



Some data were removed due to adverse meteorological conditions. 
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Table 6: Model Parameters 



Parameter 2006 Value 2008 Value 



4.2 x 10 4 3.3 x 10 4 

5.0 x 10 5 4.0 x 10 5 



Train" 
'Jmax- 

injection electron spectral index: 2.6 3.2 

Escape time parameter (t esc = rjR/c) rj esc = 3 rj esc = 3 

Magnetic field at z [G]: 0.48 0.68 

Bulk Lorentz factor: T= 20 T= 20 

Blob radius [cm] 3.0 x 10 15 3.0 x 10 15 

Observing angle [degrees] o b s = 2.87 6 Q b s = 2.87 

L e (jet) 7.76 x 10 42 erg s" 1 1.06 x 10 43 erg s" 1 

L B (jet) 3.11 x 10 42 erg s" 1 6.24 x 10 42 erg s" 1 

L B /L e 0.40 0.59 



